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■ The properties of the vacuum are addressed in the two- and four-dimensional quark 

models for QCD. It is demonstrated that the two-dimensional QCD ('t Hooft model) 



possesses only one possible vacuum state — the solution to the mass-gap equation, 
which provides spontaneous breaking of the chiral symmetry (SBCS). On the con- 
trary, the four-dimensional theory with confinement modeled by the linear potential 
supplied by the Coulomb OGE interaction, not only has the chirally-noninvariant 
ground vacuum state, but it possesses an excited vacuum replica, which also exhibits 
SBCS and can realize as a metastable intermediate state of hadronic systems. We 
discuss the influence of the latter on physical observables as well as on the possibility 
to probe the vacuum background fields in QCD. 

PACS numbers: 12.38.Aw, 12.39.Ki, 12.39.Pn 



I. INTRODUCTION 



Potential models play an important role in studies of QCD. In spite of their obvious 
shortcomings, like loss of Lorentz invariance and causality, they offer a simple and intuitive 
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tool for investigation of the chiral properties of the theory, hadronic spectra, decays and so 
on. The Hamiltonian approach, taken together with the Bogoliubov technique for diagonal- 
ization, provides a powerful method to analyse such potential models. The two-dimensional 
't Hooft model for QCD, besides being a genuine quantum filed theory, constitutes an ex- 
ample of an exactly solvable theory with an instantaneous linear potential yielded by the 
two-dimensional gluon [|TJ. 

In this paper we address the question whether QCD can possess excited vacuum states, 
hereby called replicas, on top of the usual chirally nonsymmetric vacuum. The problem 
of different vacuum states, their coexistence and stability, is not new and was addressed 
many times in many contexts (as an example, see, e.g., || or || where the "history" of 
the vacuum is briefly outlined). Here we address the theoretical possibility of existence of 
replicas in the framework of quark QCD models which have otherwise been successful in 
describing hadronic phenomenology and will argue that if one uses these sort of models - 
which indeed can be traced back to the Gaussian approximation for the cumulant expansion 
of QCD4 — for their phenomenologic success then one has also to consider the quite possible 
existence of these replicas. They come in the same package. Then it should not come as a 
surprise that the number of spatial dimensions together with the strength and form of the 
inter-quark potential should play a role of paramount importance in the existence of such 
replicas. We start from the two-dimensional QCD and demonstrate that, for two dimensions, 
only one chirally nonsymmetric state (the vacuum) may exist. Therefore in two dimensions 
there are no replicas. This statement holds true for both, linear and harmonic oscillator 
potentials. On the contrary, an infinite set of solutions exists for the mass-gap equation in 
four dimensions, if the inter-quark interaction is chosen to be quadratic in the inter-quark 
distance, as it was found in Q. In the present paper we confirm this result. However, for the 
more realistic linear potential the situation in QCD 4 changes dramatically. We study the 
corresponding mass-gap equation and find the pure linearly rising interaction to fail just to be 
sufficiently "binding" to hold any replicas, so that only one chirally nonsymmetric solution 
to the mass-gap equation may exist in this case (the trivial chirally symmetric solution is 
always present in 3+1, as opposed to the two-dimensional case). On the contrary, for the 
more physical interaction with realistic values, both for the string tension <7 and the strong 
coupling constant a s of the Coulomb interaction together with a constant term j^Uj^ , only 
to adjust for the right chiral condensate, we find one vacuum replica to exist besides the 
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usual vacuum state. We give qualitative arguments that the excited vacuum state should 
exist in the real QCD and discuss its possible influence on physical observables. 

This paper is organized as follows. In the second section we consider the two-dimensional 
't Hooft model and briefly review the formalism going through the concepts of the Hamilto- 
nian approach, the Bogoliubov transformation from bare to dressed quarks, and the mass-gap 
equation for the chiral angle. We argue that the mass-gap equation for the 't Hooft model 
has only one solution, found numerically in || , which defines the chirally noninvariant phase 
of the theory, whereas the phase with unbroken chiral symmetry possesses an infinite energy 
and, hence, is unphysical. These results are found to remain for a similar, but technically 
simpler, situation when using the harmonic oscillator instead of the linear potential. In the 
third section we deal with the case of the four- dimensional QCD and study the mass-gap 
equation numerically, both for harmonic and linear confining quark kernels. We find that 
an enumerable infinite set of replicas always exists in the case of the quadratic potential, 
whereas the linear potential appears to be just not quite enough "binding" , so that it has 
to be added either to a Coulomb or a constant potential (or to a combination of both) to 
support the existence of an excited vacuum state. It turns out that for a realistic set of 
parameters {<t , a s , U}, simultaneously describing hadronic spectra and yielding a quark 
condensate in the ball-park of (qq) = — (250MeV) 3 it was impossible to avoid the existence 
of one replica. The fourth section is devoted to the discussion of the hadronic processes in 
the presence of the excited vacuum. We give arguments on how the excited vacuum replica 
can affect the physical observables and serve as a probe of the QCD vacuum structure. Our 
conclusions and the outlook are the subject of the fifth section. Throughout this paper we 
consider the chiral limit, i.e., we always put the mass of the quark m equal to zero from the 
very beginning. 



The two-dimensional 't Hooft model, suggested in 1974 [|TJ and widely discussed in liter- 
ature as a toy model for QCD, is described by the Lagrangian density 



II. TWO-DIMENSIONAL QCD ('T HOOFT MODEL) 



A. 



Introduction to the model 




1 



(1) 
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D = (d, - igA;t a ) ltl 
and the large- A^^ limit implies that [23] 



9 2 N C 

7 = — — > const. (2) 

Air Nc^oo 

We fix the axial gauge imposing the condition Ai(xo, x) = [0, £|, so that the vector-field 
propagator takes the form 

^fao -Vo,x-y) = ~^S ab g^ g u0 \x - y|5(x - j/ ), (3) 

and integrate out gluonic degrees of freedom to arrive at the Hamiltonian, 

H = J dxq + {x) ^— il 5 ~Q~ J Q( x )~ 

Following the standard BCS technique, we introduce the dressed quark field []7|, 

q a (t,x) = J ^e^[b a (p,t)u(p) + dt(- P ,t)v(-p)}, (5) 

b a (p)\0) = d a (-p)\0) = 0, b+(p)\0) = \q), d+(-p)|0) = |g>, (6) 
{ba{p, t)b+(q, t)} = {d a {-p, t)d+(-q, t)} = 2n5(p - q)5 a p, (7) 

«(p) = r(p)Q), v(-p)=T( P )(° i y (8) 

T(p) = e-fr k >», -\<Q<\ 0{-p) = -0(p), 

where the Bogoliubov-Valatin angle 9 is subject to mass-gap equation. In the new basis the 
Hamiltonian (ffi) takes the form: 

H = LN C S V + : H 2 : + : H 4 :, (9) 

where LNc£ v stands for the vacuum energy (L being the one-dimensional volume) and the 
second and the third terms on the r.h.s. for the quadratic and quartic in the quark/ antiquark 
creation/annihilation operators, respectively. In the remainder of this paper we shall con- 
centrate on the energy density S v . 
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B. The vacuum energy and the chiral symmetry breaking 

The vacuum energy LNq£ v can be calculated as an average of the time ordered Hamil- 
tonian over the BCS vacuum state [24|. 



E vac = (0\TH(b, b + , d, rf + )|0) = LN c £ vac = LN C (£%? + A£ vac ), (10) 

where we have introduced the excess of the vacuum energy density over the free-theory 
energy density, £J^?. We have 

A£ vac = / ^Tr | 7 5 pA-(p) + ^ / J^ A +( k ) A -(P)} ~ ("J 
A±(p)=T(p)^pV(p) 

or, in terms of the angle 9, 

M™[0] = ~J ^(P^9(p) -\P\)-^J j^y 2 cos [B(p) - 9{k)) , (12) 

where for the free massless theory we substituted 7 = and 9f ree (p) = |sign(p). 

The excess A£ vac [#] is the main object of our study. The true angle 9 should minimize 
the vacuum energy, that is, it is a solution to the mass-gap equation, 

First, we perform a simple qualitative analysis of Eq. (|I2"D . Let 9(p) be the solution of ( [T~3[ ) 
corresponding to the minimum of A£ vac [9}. Then for 9(p/A), with an arbitrary parameter 
A 1, A£ v3iC [9/A] > A£ V3JC [9]. Naive dimensional analysis shows that A£ vac scales with A 
as 

A£ vac = ^jCiA 2 - 7C2 (14) 

with C\p being positive constants. Then the stable solution is provided by minimizing the 
energy fll4p with respect to A, that is, by A = 0, which corresponds either to the free 
massless theory, 

7T 

0free(p) = -sign(p), ^freeb) = \p\, (15) 

with Ef iee (p) being the free-quark dispersive law, or to the analytical solution of the mass-gap 
equation found in [^], which reads: 

9(p) = ^sign(p), E(p) = \p\-P^- (16) 
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where the symbol P stands for the principal-value prescription. Both solutions yield un- 
broken chiral symmetry. Thus one would have arrived at the erroneous conclusion that no 
chirally nonsymmetric solution of the mass-gap equation may exist, were not for the fact 
that the vacuum energy (O) is logarithmically infrared divergent and the correct form of 



the relation (|TJ|) should have been 

A£ vac = X -C X A 2 - 7 C 2 In A + 7 C 3 (17) 

instead of Eq. (|bj|) , with C3 proportional to the logarithm of the cut-off p5[| . The function 
(|T7|) always has a minimum at 

A = ^ 0, (18) 

which corresponds to a nontrivial solution of the mass-gap equation, found numerically in || . 
From the form of the function ([lj]) one can see the logarithmic growth of the energy when 



approaching the solution (|T6|), which corresponds to A = p6 |. Therefore two conclusions 
can be deduced. First, there is no need to demand A£ vac to be negative as the nontrivial 
vacuum energy defined by (0) is always compared with the infinite energy of the chirally 
symmetric phase. Second, no phase transition of chiral symmetry restoration is possible in 
the 't Hooft model, its phase diagram being trivial. A similar conclusion is made in [|1(J . 
The chiral condensate (qq) is given by 

N c 



7T 



dp cos 6 (p). (19) 



Then the nontrivial solution, found in |J, provides a nonvanishing chiral condensate, whereas 
the solution ( |T6| ) gives (qq) = 0, but with an infinite energy. Thus we conclude that in the 
't Hooft model there can only be one phase, with spontaneously broken chiral symmetry. 
Note that there is no contradiction to the Coleman theorem ||TTJ, which precludes sponta- 
neous breaking of symmetries in two-dimensional theories, if one considers infinite- Nc limit 
I PI- 



The contribution of the infrared logarithmic divergence in (|17|) has yet another conse- 
quence for the 't Hooft model — namely, that the nonvanishing chiral condensate and the 
spontaneous chiral symmetry breaking happen already at the BCS level, whereas going be- 
yond BCS just provides corrections suppressed by powers of Nc- Indeed, staying at the BCS 
level, one diagonalizes the theory in the one-particle sector, introducing a nontrivial vacuum, 
full of correlated qq pairs. It was demonstrated in M that one can go beyond BCS in the 
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't Hooft model performing the second Bogoliubov-like transformation and diagonalizing the 
theory in the mesonic sector, that is, in the sector of the quark- ant iquark bound states. The 
true vacuum of the theory, is connected to the BCS vacuum, |0), via a unitary operator. 
One should use the true vacuum state for calculations of all matrix elements of operators, 
like the chiral condensate, but as shown in ||, the difference between the two averages is of 
the next-to- leading order in Nc' 

(Q\qq\Q) = (0\qq\0) + O (v^) ~ 0{N C ) + O (v 7 ^) • (20) 
It is instructive to compare the 't Hooft model with other two-dimensional models, e.g., 



with the Gross- Neveu one |[L3|| . The latter is given by the interaction Lagrangian 

1-2 

Lint = -j9 (4>i4>i) , 9 > 0, (21) 

where the flavour index i runs from 1 to > 1, whereas A = g 2 N remains finite. The 
positive sign of the constant g chosen in fl21"|) is known to lead to condensation in the singlet 



ipip channel Jl3| , whereas the negative sign leads to charged ipip and ipip condensates [14J . It 
was shown that the chiral symmetry is also broken in the Gross-Neveu model, but in order 
to have a nontrivial chiral condensate one has to go beyond BCS level, summing fermionic 
bubbles. Bosonization of the model, in terms of the compound state a = gipip, appears to 
be the most convenient formalism for doing this and the resulting renormalized "potential" 
reads: 

1 \ F / rr \ 2 

(22) 



ln[ *Y 



with er being the renormalization point, which brings the mass scale into the theory. Only 
the first term in ( p2|) appears at the BCS level leading to the trivial solution <7 m i n = 0. On 
the other hand, the full potential (E21) provides also a nontrivial solution, 



cr min = ±cr exp (l-^j . (23) 

From Eq . (^2|) one can see that the logarithmic term, which scales as In A under the above 
mentioned transformation, is multiplied by o 2 ~ A 2 which cancels the logarithmic growth 
of the potential at the origin. Thus the two phases, chirally symmetric and nonsymmetric, 
coexist in the model, whereas the latter, with nonvanishing <7 m i n , is energetically preferable. 

The chiral condensate of the Gross-Neveu model scales like 

(a) ~ (gipip) ~ gN ~ \\/N ~ \/N (24) 
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for large N. Relation ( p4[) is analogue of (|20"D with vanishing leading term of order O(N), 
which is in agreement with the statement made above, that the chiral condensate in the 
Gross-Neveu model vanishes at the BCS level. One can arrive at the same conclusion from 
simple dimensional considerations. Indeed, the Gross-Neveu model describes scale-invariant 
renormalizable theory, so that there is no parameter with dimension of mass in the bare 
Lagrangian and one has to deal with the mechanism of the dimensional transmutation to 
get a scale. In contrast, the 't Hooft Lagrangian (|l|) defines a super- renormalizable theory 
and the dimensional parameter, the coupling constant g, is present in the theory from the 
very beginning. 

Formula can be given a more physically transparent form if one notices that the 
chiral condensate X = (qq) from Eq . (|19D scales linearly with A, X — > AS, so that it can be 
used instead of A, 
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3; 



(25) 



where S n stands for the actual value of the chiral condensate. 



C. Harmonic oscillator-type interaction 

In this subsection we complete the study of two dimensional Hamiltonians using a har- 
monic oscillator-kernel in @, instead of the linear kernel, 



or, in momentum space, 



1 \x-y\ ^K^(x-y) 2 , (26) 



7 _^3r« 



{ji _ h)2 .«K*S"(p-k), (27) 

with Kq being the new constant with the dimension of mass. Due to the presence of the 
delta function in k in the interaction kernel, all integrals can be considerably simplified. For 
example, the excess of the vacuum energy A£ vac [6>], similar to (0), reads now: 

AS vac [6] = - J g |[psin%) - b|] - ^ [6'{p)f} . (28) 

The integral on the r.h.s. in (p8|) is convergent and scales as 

A£ vac [#] = C,A 2 + (29) 
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FIG. 1: The ground-state solutions of the two-dimensional mass-gap equations ( p2| ) which corre- 
spond to the linear (left plot) and the oscillator (right plot) potential, respectively. The momentum 
p is given in units of yfiy for the left plot and in units of Kq for the right one. 



under 9(p) — > 9(p/A), so that qualitatively the same conclusion holds true, that there exists 
only one phase of the theory with spontaneously broken chiral symmetry and there is no 
unbroken phase. Similarly to (p5|), we can rewrite Eq.(|29|) using the chiral condensate instead 
of A: 

2 V£r 



(30) 



with the minimum given by E = E 



o- 



D. Numerical solutions for the chiral angle 

In the previous subsections we found that there must be nontrivial functions 9(p) pro- 
viding the minimum of the vacuum energies ([12]) and (p8|) . As was mentioned above, these 
functions are solutions to the mass-gap equations, 

SAS mc [9] _ 5AS^9]_ 

69(p) ~ U ' 89{p) ~ U ' {6l) 

which take the following forms: 

pcos9(p) = lj ^_sin[9(p) - 9(k)), pcos9(p) = -^9"{p) (32) 
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FIG. 2: The unphysical excited solutions to the mass-gap equation ( |32| ) (oscillator confinement). 
The momentum p is given in units of Kq. 

for the linear and harmonic-oscillator potential, respectively. 

Now we turn to numerical studies of these two equations. Note that the qualitative 
analysis performed in the previous subsection cannot answer the question of how many 
solutions of the mass-gap equation exist. Indeed, a second solution, if it exists, would just 
lead to different constants C"s (C's) and a different value of the chiral condensate Eo (So), 
since this solution is expected to have a different profile and it cannot be obtained from the 
ground-state solution by dilatations. 

As the first step, we reproduce the solution of the first equation in (|32| ) found in and 
construct a similar solution for the other (see Fig. 1). Notice that we define the chiral angle 
9(p) from — 7r/2 to tt/2 with the boundary conditions 9(p = 0) = 0, 9(p — ► oo) — > tt/2 (as 
was mentioned before, 9(—p) = —9{p)). It turns out that under such conditions no excited 
solutions of the mass-gap equations fl32[) exist. Only if both requirements are relaxed, can 
one find the solutions depicted in Fig. 2, which start at for p = and tend to tt/2 + 27m, 
n > 1, for p — > oo. From Fig. 1 one can easily see that the solutions for the ground state 
in both types of the confining potential have a very similar shape, so it would be natural 
to expect similar excited solutions to exist for the linear potential as well. However, such 
solutions lead to maximal condensation of the qq pairs not only at p — 0, but for some finite 
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p's, where cos 6 = 1, which is unphysical, so we disregard these solutions. 

Thus we conclude that only the ground-state solution may exist for the mass-gap equation 
in the 't Hooft model, defining the only possible phase of the theory, with spontaneously 
broken chiral symmetry and the chiral condensate 

S = -0.410iV Cv ^, (33) 

which coincides with the value found in a sequence of papers || § ||. A similar result for 
the theory with the oscillator-type interaction reads: 

S = -0A05N C K , (34) 

which is very close to ( |33"D if the scales y/j and K are chosen close to one another. 

Notice that in the two-dimensional QCD the energy of the vacuum appears to be an 
apparently meaningless quantity by itself, although it leads to the meaningful mass-gap 



equation, as discussed above. Indeed, the expression (|TJ) contains a divergent integral, so 
that the result depends on the cut-off. On the other hand, in view of the fact that the 
physical solution of the mass-gap equations is unique, one cannot define a difference of two 
energies, where the divergent term would cancel. 

Now, when the numerical methods are tuned and tested at the simple example of the two- 
dimensional QCD, we turn to the main subject of the present paper — namely, to studies 
of the vacuum replicas in QCD 4 . 



III. FOUR-DIMENSIONAL QCD 



A. Instantaneous interaction and forms of the potential. The mass-gap equation 

In this subsection we consider the case of QCD 4 and study the problem of existence of 
the vacuum replicas for various types of gluonic correlators. In contrast to QCD 2 , we put 
the number of colours equal to three from the very beginning. 

Following the standard procedure, we start from the QCD partition function and perform 
its cluster expansion, so that, after integrating out the gluonic field, the theory contains 
effective 4-quark, 6-quark, and so on vertices with the formfactors given by the corresponding 
correlators of the gluonic fields, 



12 



where the irreducible average, ((...)), is defined in the standard way, 

«1» = <1>, «12» = (12>-<1><2>, .... (36) 

We fix the Fock-Schwinger gauge (x — ZQ)^A a {x) = with an arbitrary fixed point z 
that allows us to connect the gluonic field with the field-strength tensor, F®(x), 

Afa) = [ ada{x - z ) u F^{a{x - z ) + z ), (37) 



and, finally, to deal with the set of gauge-invariant cumulants [fL7|, 



(^3,^0))), 



(38) 



F» (x,Zq) = $(zo,x)F*(x)$(x,Zo), 



with $(x,y) being the standard parallel transporter from the point x to the point y. Now 
we apply a set of approximations. 

On one hand, the Gaussian dominance is used, which means that we keep only the low- 
est-order correlator ((FF)), whereas all others are considered suppressed []I7| . Such an 



approximation can be justified, first, by the Casimir scaling, suggested long ago [IS] and 



recently confirmed by lattice calculations |1| . The Casimir scaling is exact in the Gaussian 



approximation. On the other hand, the minimal-area law asymptotic for an isolated Wilson 
loop can be saturated by the Gaussian correlator. Finally, numerical calculations of hadronic 
observables within this method show good agreement with experimental and lattice data. 
The details of the method can be found, e.g., in the review 0] and references therein. 

As far as the form of the nonperturbative part of the bilocal correlator is concerned, we 
choose it in the form 

((A a »(x)Al(y)}) 
1 ,1 

ada / (3d(3(x - z ) x (y - z ) a ((F^(a(x - z ) + z )F^((3(y - z ) + z ))) (39) 
Jo 

-> S^g^gMdx - y\)5(x - y ), 

which implies that the gluonic correlation length T g (also playing the role of the QCD 
string radius) is small enough, so that the string-like interaction can be approximated by 
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the instantaneous one [p7| . For the confining potential, Vo(|x|), we consider two cases, the 
quadratic and linear form, 



V Q {\x\) = Kl\x\\ 
V (\x\) = a \x\ 



(40) 
(41) 



The Lorentz structure of confinement is taken to be 70 x 70 . In the linear case we shall add 
to Vb(|a?|) a Coulomb term and a constant, 



a. 



V 1 (\x\) = --^ + U, 



x 



(42) 



with the Lorentz structure being 7^ x 7^. Notice that due to the colour structure of the 
interaction, being ^- x the constant term in the potential fll"2]) yields a force between 
singlets which is proportional to the Resonating Group Method normalization overlap of 
these singlets. This is a simple example of the Wigner-Eckart theorem. The Coulomb 
potential is necessary for this class of models to possess a heavy-quark limit compatible 
with the known spectroscopy. 

Thus our model is defined by the parameter set {o"o, a s , U}. We shall choose the usual 
literature values for two of them, o~o = 0.135Ge\^ 2 , a s = 0.3, and adjust U to obtain a 
quark condensate (qq) ~ — (250MeV) 3 . Next we turn to the study of the vacuum structure 
following the same steps used in QCD 2 . 



First, we introduce the dressed quark field |L5], [22 

q a {t,x) 



d 3 p 



E J ^e^[6 aC (p,t) Wc (p) + d+(-p,0^(-P)]> 



(43) 



with a and ( being the quark colour and the projection of the spin, respectively, and the 
amplitudes u and v defined as 



u(p) = 
v(-p) 



l+sin tp 



+ 



1— sin</3^ -> 

— Y^P ■ Oi 



U ( 



l+sin ip 
2 



(44) 



Here the chiral angle ip is used instead of QCD 2 9. They are simply related by 

<p(p) = l-9(p). (45) 



It is convenient to define the new angle, if, from — 7r/2 to 7r/2 with the boundary conditions 
f(p = 0) = 7r/2 and cp(p — > 00) — > 0. 
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If the chiral angle <p(p) is not identically zero, then the true vacuum state is populated 
by the quark-antiquark pairs with the 3 Po coupling, 

|6) = n dp ( 1 + C0 2 S2 ^ + ^CJ p + ^fM C g) |0) . S |0), (46) 

where indices c and / numerate the colour and the flavour, respectively; the angle 4>{k) is 
connected to the chiral one, (p(k), via a simple relation, 



m 



(pip) = arctan h 2(f)(p) = 2(p(p), 

p m=0 

and defines the above mentioned unitary operator S . 



(47) 



C+ p =(b+(p),b+(p)) cf dJl 



d+(p) 



(48) 



^, = (-V6)E 1 10 k 1/2 1/21 H^- A ' 

CTjm y m a o y y en o- 2 o- y y p 2 p x + w y 

The operators annihilating the new vacuum, |0), are 



(49) 



Kp) = S b(p)S£, d(-p) = S d(-p)S£. 



(50) 



The interested reader can find the details of the formalism in papers p2| . 

The next step consists of arranging the normal ordering of the QCD 4 Hamiltonian, 



H = I d' 6 x[H {x) + Hf\x) + HY\x)} 



r( 2 )/ 



(51) 



H (x) = q + (x) (-ia ■ V J q(x) 



(52) 

Hf\x) = \ I d 3 yq + (x)^-q(x) [v (\x - y\) + V 1 {\x - y\)] q+(y)^-q(y), (53) 



Hf\x) 



1 C X a \ a 

d 3 y q + {x)d—q{x) V x {\x - y\) q + {y)a—q{y). 



(54) 



2 J " " y ' 2 " v ' Vl °" " 2 
In the new basis, we obtain an expression similar to @ and we again concentrate on 
the minimization of the vacuum-energy density — namely, on the mass-gap equation, which 
takes the form: 

A{p) cos ip(p) — B(p) smip(p) = 0, (55) 



A{p) 



d 3 k 
(2^)3 



V (p - k) + 4Vi(p - k) sinip(k) 
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S(p) = p + -C F J —— 3 {p- k) [V (p - k) + 2T4(p - k)\ cos 

where Cf = § is the Casimir operator in the fundamental representation. 

Thus, equipped with the general form of the mass-gap equation, we are in the position to 
perform the numerical search for the excited vacua for the oscillator and linear potentials. 



B. Solutions of the mass-gap equation for the oscillator-type interaction 

In this subsection we solve numerically the mass-gap equation fl55"D for the oscillator-type 
confining potential (f40D and V\ put to zero. Similarly to the two-dimensional case, the 
Fourier transform of the potential (fl0| ) is the second derivative of the delta function, 

Vo® = -WK 3 ^\p), (56) 

which leads to the mass-gap equation in the differential form, 

2 2p 1 

V" = — V?'+ n rs 3 siny- — sin2^, (57) 
p C fKq p z 

and we solve it numerically using the Runge-Kutta method. In Fig. 3 we show the ground- 
state and the first excited solution of Eq. (|57|) . 

Thus we conclude for the existence of at least one replica for the case of quadratic confining 
potential. Moreover, using the same numerical technique, one can build the third, the fourth, 
and so on excited solutions, as was found in 0]. In fact, the harmonic oscillator supports 
an infinite tower of such replicas. 

C. Vacuum replica for the linear potential 

Finally, we turn to the main goal of the present paper — namely, whether the linear 
confining potential supports the existence of replicas. 

In contrast to the case of the oscillator potential, the mass-gap equation for the linear 
confinement, given by Eq.(il), is integral, 

psiny?(p) = ^ ° f ^— sin Lf(k) cos tp{p) — p- k sin ip{p) cos (f(k) , (58) 

27T J (j) — k) 4 L -I 
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FIG. 3: The first two solutions of Eq.(|57|) for the chiral angle <p(p) found with the Runge-Kutta 
method. The momentum p is given in units of Kq. 



or, equivalent ly, 

u f°° 

psinipOp) = — / dk k z 
4tt J 

f2k(k 2 +p 2 



(p 2 — k 2 ) 



- sin cosip{p) 



(59) 



In 



p + k 



cosip(k) sin ip(p) 



\p{p 2 — k 2 ) 2 \r 1 1 
where a = Cf&o is the string tension in the fundamental representation, and the angu- 
lar variables are integrated out. One can check that the integral on the r.h.s. of (|59D is 
convergent. Eq.(p9|) possesses only one nontrivial eigenstate: the ground-state vacuum. Ev- 
erything happens as if the confining linearly rising potential just fails to be "binding" enough 
to produce an excited vacuum replica. Besides, the chiral condensate appears strongly un- 
derestimated in case of the pure linear confining potential . Then we supply it with the 
extra interaction, given by (f42D, which reads, in momentum space, 

47TOU 



Vi(p) 



P 2 



(2ir) 3 U6®(p), 



and modify the mass-gap equation accordingly, 

4 



p sin (p(p) 



dk k z 



a 

/ 2(k 2 + p 2 



(p 2 - k 2 ) 2 
) 1 



sin tp{k) cosy?(p) 
p + k 



(60) 



(61) 



\pk(p 2 — k 2 ) 2 p 2 k 2 



In 



p — k 



cos (f(k) sm(f(p) 
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CpOl s 

4:71 



{ P 2 + k\ n 

p 2 k 2 



/■OO 

/ dk k 2 




p + k 


/o 


_pk 


p — k 



p + k 



P 



smip(k) cos ip(p) 



2 \ 

— cosip(k) smipip) 
pk J 



+-C F U sin 2tp(p). 



In contrast to Eq.(^), the mass-gap equation (^T|) is not divergence- free because of the 
Coulomb part, proportional to a s . If we considered the case of a nonzero quark mass m, 



then the given divergence would be only logarithmic, however, for Eq.(61) we have m = 
and the divergence becomes linear. We regularize it by using a modified potential Vi, 

1 1 



Vi(p) = 4ira s 



p 2 



+ (2tt)W 3) (p), 



(62) 



_p 2 + A 2 

with a cut-off A. Notice that the mass-gap equation for the single Coulomb potential does 
not contain any scale (it would be given by the quark mass m beyond the chiral limit), 
which appears only after the regularization and is given by A. Then all quantities with the 
dimension of mass are directly proportional to A. 

Now let us assume, for a moment, the general case of an arbitrary nonvanishing current 
quark mass. Then the quarks acquire an effective mass m e g, which can be written as a sum 
of two components, 

m c fj = m + E(m, A), (63) 

where m and E(m, A) are the current quark mass and the quark selfenergy due to the 
contribution of the chiral condensate, respectively. The latter is obviously A-dependent. To 
evaluate m eff one can use the relations 

m(p) 



sin ip(jp) 



cos ip(p) 



P 



(64) 



\/p 2 + m 2 (p) ' \fp 2 + m 2 (p) ' 

similar to the free-case relations, but with (pip) being the solution to the mass-gap equation. 
Now the effective mass, m e g, can be associated with m(p ~ 0). Let us consider the heavy- 
quark limit, Tn — ► oo, first. A large quark current mass is known to destroy the chiral 
condensate. As a result, the self-energy part of the relation (B3) must vanish, so that the 
effective mass is basically given by the current mass, 



m e fj ~ m. 



(65) 
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TABLE I: Parameters of the model fixed from the fit for the ground vacuum state. 



a , GeV 2 


a s 


U, MeV 


A, MeV 


m e fj, MeV 


(qq), MeV 3 


0.135 


0.3 


220 


250 


230 


-(250) 3 



In particular, £(m, A)/m — > 0. That is, for large m, we have to go up in A for the system 
to become more and more Coulombic, whereas the chiral condensate goes to zero. In other 
words the potential of Eq.(|62|) naturally possesses the Coulomb limit for very heavy current 
quark masses. 

In the opposite limit, m = 0, the effective mass appears entirely due to the chiral sym- 
metry breaking, 

m eff = E(0, A) = A x Const, (66) 

where the last equality holds true provided we choose A ~ y/&o, the only scale of the 
problem. Thus in either case m e g > A, and it is remarkable that such a potential should 
equally address both heavy- and light-quark sectors. 

Therefore we adopt the following strategy. We choose the standard values of the string 
tension, a = 0.135GeV 2 , and the strong coupling constant, a s = 0.3. As discussed above, 
we put A to be of order of ^/Wq and solve the resulting mass-gap equation numerically for 
the ground state, fitting the only free parameter, U, in such a way, as to have the chiral 
condensate, 

3 z* 00 

(qq) = 2 dpp 2 sinp(p), (67) 

71 Jo 

to be around its standard value, about — (250MeV) 3 . With the solution found for the chiral 
angle tp(p), we evaluate the effective quark mass, m e s, and check if it has the value of order 
200 -T- 300MeV, that is, the value of the same order of magnitude as the regulator A, being 
also close to the value usually adopted for the constituent quark mass. For selfconsistency 
we also need to ensure that the fitting parameter, U, also appears to be of order of the 
interaction scale and, hence, does not bring a new scale into the problem. In Table 1, besides 
the known values of do and a s , we give the best set of parameters U and A consistent with 
the standard value of the chiral condensate. 

Now, when all parameters of the model are totally fixed, we perform our search for the 
excited vacuum replica which, indeed, appears to exist, and we give its profile, together 
with the ground-state one, in Fig. 4. We tried many other combinations of U and A and, 
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2.5 



FIG. 4: The ground-state (left plot) and the first excited (right plot) solution of Eq.(|6l|) for the 
set of parameters given in Table I. The momentum p is given in units of yTro. 

whatever values we have used, we always found this replica to persist. 
The chiral condensate for the excited vacuum state, 

(qq) = (12QMeV) 3 , (68) 

is smaller in the absolute value and has a "wrong" sign compared to the ground-state 
one, given in Table I. For massless quarks this sign can be easily reversed by the change 
ip — ► — ip everywhere, which also satisfies the mass-gap equation and does not affect the 
vacuum energy. However, nonvanishing quark mass breaks this symmetry, since an extra 
term, mcos(p(p), appears in the mass-gap equation ([58]). It follows immediately from the 
Gell-Mann-Oakes-Renner relation that the sign of the condensate given in (|68| ) leads to the 
pion with an imaginary mass, that is, to the tachyon. However, the latter does not lead to 
a disaster since the excited vacuum state is, indeed, unstable and it tends to decay to the 
ground state. 

Now we calculate the density of the vacuum energy defined as 

A^vac = -f / 7^8 { A ^ Sil1 V(P) + W + Pi C ° S V(P)) - £ - " ( 69 ) 

where A(p) and B(p) are given in (|55|) and the degeneracy factor g counts the number of 



20 



independent quark degrees of freedom, 

g = (2s + l)N c N f , (70) 

with s = | being the quark spin; the number of colours, Nq, is put to three, and the number 
of light flavours, Nf, is two. Thus we find that g = 12. The trivial solution of the mass-gap 
equation is ^trivial (p) = (in contrast to the two-dimensional case, this solution does not 
lead to difficulties with the vacuum energy). Then for the two nontrivial solutions, depicted 
in Fig. 4, one has 

ground = _( 128 Me\/) 4 (71) 

and 

A £excited = _( 36Me ^)4. ( 72 ) 

The negative sign ensures that they both are more energetically preferable than the trivial 
vacuum. 

It is instructive to rewrite Eqs.([H]), ( |72]) in the form of the vacuum energy per Fermi 
cubed: 

A £gTound = _ 3 4.7Me^//m 3 , (73) 
^excited = _ .2MeV//m 3 . (74) 

Now, if the vacuum is locally excited, then the extra energy stored in the volume Vo = 
lfm 3 is 

AE = (A£™f cd - A£ v s a r c ound ) Vo « U.hMeV. (75) 

We have checked the result ( [75] ) against slight variations of the parameters of the model 
and found it to be very stable, changing only by several per cent. 

IV. PHYSICAL PROCESSES IN PRESENCE OF THE EXCITED VACUUM 

REPLICA 

In this section we show how the excited vacuum could be "seen" in hadronic processes. 

As follows from Eq. ([75]) , the energy stored in the bubble with the excited vacuum inside 
is proportional to its volume. At the moment we can only give hand-waving qualitative 
arguments concerning the most probable value of this volume. Indeed, the energy scale of 
the interaction we have in QCD is about 300MeV that corresponds to the distances about 
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0.5/m which can be expanded up to the Compton length of the pion, being about 1.5/m 
or so — the standard hadronic scale. With the bubble of such a radius one can achieve 
the energy around 300MeV, which is sufficient to produce an S-wave pair of pions when 
the excited vacuum decays back to the ground-state one. This value looks quite naturally 
since all parameters of the model, like the interaction scale, y/&o, and the SBCS one, given 
by the chiral condensate or by m e g (see Table I) have the same order of magnitude. As 
was discussed above, all three scales should always appear selfconsistently, regardless of the 
details of the model used for numerical calculations. In the present paper we give just an 
example of such a model, which we find rather realistic. 

It is instructive to see how the excited vacuum enters the amplitudes of the hadronic 
processes. Usually one is to evaluate the matrix element of the form 

(0, hadrons\O\0, hadrons') (76) 

with O being a field operator. Note that this is the true vacuum of the theory, |0), to 
enter the matrix element, and we put it explicitly in (|76|) . If the excited vacuum comes 
into the game, then it can appear either in the bra or in the ket vector, or in both. We 
use the notation |1) for it. Let us illustrate this statement by the correlator of two quark 
electromagnetic currents, 

(0\TJ^x)J u (y)\l) = (0\Te ip *J^{0)e- ip *e ip yU0)e- ip y\l), (77) 

where P is the operator of the total momentum. Now, disentangling the T product and 
making the external operators e ±iPx and e ±lPy act on the vacuum states to the left or to the 
right, one can arrive at the expression 

[ e- iAPy (0\ JJx - y) J„(0)|1), x > y 
[ e-^iOlJ^y - x) J„(0)|1), y > x , 

where AP is the difference of the four-momenta of the two vacua, which, actually, has only 
the zeroth component — the difference of the vacuum energies, AE. If the Fourier transform 
is now applied to (|75|), 

V^(p,q) = j Mye^e^(tyTJ^x)JM$), ( 79 ) 



then 



V^(p, q) ~ 5(po - go - AE)6^ 3 \p - q), (80) 
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that is, one has an effective source pumping the energy into the system. This energy can 
be spent to produce photons on mass shell each with an energy u = AE/2 ~ lOMeV for 
the radius of the excited- vacuum bubble being of order of the interaction radius, that is, of 
order about 0.5/m. 

Notice that if the BCS vacuum structure is ignored and the trivial vacuum, |0), is used 
instead of the |0) and |1) in ( |77| ) (this corresponds to the chiral angles ip and (p identically 
put to zero), then the latter expression coincides with the standard perturbative- vacuum 
polarization operator supplied by the external energy source. In order to perform calcula- 
tions taking into account the nontrivial functions ip and (p, one should express all quantities 
in terms of the quark creation and annihilation operators inherent to one and the same vac- 
uum, and to apply the standard Wick contractions then. For example, relating everything 
to the trivial-vacuum operators, one can find: 

|0)=£„|0), |l) = Si|0) (81) 

with operator So defined in fllSD and a similar definition of S\ with the obvious change 
ip — > (p. Then the correlator ( |78"|) reads: 

(0\TJ^)My)\l) = <0|S; ^:^f Sl|0> . (82) 

where we used the correct normalization which removes the disconnected diagrams. Note 
that in flS2|) the currents are also expected to be expressed in terms of the trivial-vacuum 
quark creation/annihilation operators, so that, starting from the quark field (|43|), one must 
use transformations fl50f) (and similar ones for the excited vacuum replica) to have everything 
selfconsistently. 

An important comment is in order here. In the above analysis we did not consider the 
mechanism of the vacuum excitation from the ground state to the replica. Indeed, to excite 
the vacuum inside of some region in space, one should have a trigger — a localized operator 
with the vacuum quantum numbers. Such an operator could be constructed out of the 
background gluonic fields (see Eq.fl38|) above) - the only bricks at our disposal. If so, the 
outgoing photons or pions, produced as a result of the excited vacuum decay, have to carry 
information about the vacuum-gluonic-field correlations and, if measured experimentally, 
offer a possibility of the direct probes of the QCD vacuum structure! 
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V. CONCLUSIONS AND OUTLOOK 

In the present paper we study the properties of the vacuum states in the potential models 
for QCD. We find that there is only one vacuum in the two-dimensional case, and it provides 
spontaneous breaking of the chiral symmetry. Lack of extra dimensions prevents this simple 
model from having vacuum replicas. On the contrary, the four-dimensional theory reveals a 
richer vacuum structure, not only containing the trivial solution with unbroken chiral sym- 
metry, but also possessing excited vacuum states, which we find numerically for the oscillator 
and linear confining force. In the latter case the Coulomb and the constant potential were 
also added to have a better correspondence with real QCD. Thus we draw the following 
picture. Starting from the theory with bare quarks and performing the Bogoliubov-Valatin 
transformation to "dress" them, we find several different types of this "dressing", which 
correspond to different vacuum states and differently broken chiral symmetry. The mostly 
broken phase, that is, the one with the maximal (in the absolute value) chiral condensate and 
the solution for the chiral (Bogoliubov) angle without knots, possesses the minimal vacuum 
energy and, thus, defines the true vacuum of the theory. Besides, there are other solutions 
of the mass-gap equation with one, two, and more knots, defining other phases of the theory, 
which are less energetically preferable than the ground state, but still possess lower vacuum 
energy than the trivial vacuum. All these excited solutions lead to spontaneous breaking 
of the chiral symmetry, though the chiral condensate decreases (in the absolute value) for 
each next replica of the vacuum. The number of such replicas depends on the parameters 
of the interaction, and we found only one of them to exist for the combination of the lin- 
early rising, Coulomb, and constant potential with the realistic set of parameters, listed in 
Table I. We tried several contributions of the the parameters a , a s and U consistent with 
the phenomenology and numerically found it to be impossible to get rid of this replica. 

Thus starting from the system in its true vacuum state and exciting it, e.g., by means of 
heating or placing a source of the strong field, one will finally arrive at the totally disordered 
chirally symmetric phase, though at the intermediate stage the system will quasistabilize 
(one or several times) that corresponds to its reordering over a new, excited, vacuum state. 
It is no surprise that all excited vacuum replicas are metastable since the hadronic spectra 
build over them contains tachyons. Still they may live for a finite period of time comparable 
with the characteristic one of hadronic processes. Then the latter may go through the 
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intermediate stage with formation of the local bubble of the excited vacuum which decays 
then to the ground-state true vacuum emitting pairs of photons and/or pions. The latter 
could be detected experimentally and, thus, serve as a signal of the existence of the excited 
vacuum state. Moreover, such measurements can serve as a direct probe of the QCD vacuum. 

Of course, the question may be raised, to what extent we rely on potential models of QCD, 
and if the statement concerning the excited vacuum replicas is, indeed, model independent. 
However potential models were shown to be very efficient in studies of chiral symmetry 
breaking, hadronic properties and decays, and so on. Thus we do believe that, at least 
qualitatively, their predictions, including the existence of excited vacuum replica, will persist 
in more realistic models for QCD, as well as in the true theory, if it is ever solved. 
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